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ABSTRACT 



Aims. The aim of this paper is to investigate how the new imaging Polarimeter and Coronograph (PolCor) at the Nordic Optical 
Telescope* (NOT) can be used in the study of circumstellar structures around AGB stars. The purpose is to prepare for a study of a 
larger sample. 

Methods. We have observed two types of AGB stars using the PolCor instrument on the NOT: the binary S-type star W Aql and 
two carbon stars with detached shells, U Cam and DR Sen The polarized light traces the dust distribution around the stars. From the 
polarimeter images the polarized intensity, the polarization degree, and the polarization angle over the images are calculated. The 
location and extent of dust structures are examined in the images. The total dust mass and the dust-to-gas ratios of the detached shells 
are also calculated. 

Results. The images of the circumstellar envelope of W Aql show what seems to be an elongated structure in the south-west direction. 
The detached shells of U Cam and DR Ser are clearly seen in the images. This is the first time the detached shell around DR Ser has 
been imaged. The radii (Rsh) and widths (AR,!,) of the shells are determined and found to be R^\^=l'!9 and 7'.'6, and A/?sh=0'.'9 and V!2, 
for U Cam and DR Ser, respectively. This is consistent with previous results. The dust masses of the feature south-west of W Aql, and 
in the shells of U Cam and DR Ser are also estimated and found to be 1 x 10"*, 5 x 10"^, and 2 x 10"* Mq, respectively. 
Conclusions. W Aql is a known binary and the shape of the circumstellar envelope seems to be in line with what could be expected 
from binary interaction on these scales. For the shells, the results are in agreement with previous investigations. Ages and formation 
time-scales are also estimated for the detached shells and found to be consistent with the thermal-pulse-formation scenario. 

Key words. Stars: AGB and post-AGB - Stars: imaging - (Stars:) binaries: general - Stars: carbon - Stars: mass-loss - Stars: late-type 
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1. Introduction 

All stars with masses between -0.8 and SM© will eventually 
evolve up the Asymptotic Giant Branch (AGB). The life time on 
the AGB, the nucleosynthesis, and the amount of dust and gas re- 
turned to the interstellar chemical cycle are all strongly affected 
by the mass-loss rate of the star during this phase. This makes 
the mass loss the most important process for the final evolution 
of low- to intermediate-mass stars (Herwig 2005]). 

The processes governing the mass loss of the AGB stars 
are not understood. In general, the mass loss is assumed to 
be smooth, spherically symmetric and driven by a pulsation- 
enhanced dust-driven wind. However, recent advances have re- 
vealed what appears to be a more complicated picture. 

Images of light scattered by the circumstellar dust have 
revealed arcs, elongated and bipolar structures, and even 
spiral shapes around a number of well-known AGB stars 
dMauron & HugginsI l2000l |2006|) . Whether this is a result 
of how the matter is expelled from the s tar or late r inter - 
action in the wind is under debate (e.g., iLeao et al] l2006h . 
Reports on observations indicating a clumpy circumstellar 
medium (both gas and dust), become more common as the 



resolution of the observations increase (IWeigelt et all l200l 
Schoier etaH l2006t ICastro-Carrizo et al.l l2010riOlofsson et all 
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* Based on observations made with the Nordic Optical Telescope, op- 
erated on the island of La Palma jointly by Denmark, Finland, Iceland, 
Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Institute de Astrofisica de Canarias. 



2010i)- Interferometric observations of OH, SiO and H2O 
maser emission also indicate a clumpy gas distribution and, 
in some cases, w hat appears to^ be bi polar outflows or jets 
close to the stars (Diamond et al. '1994': 'Szvmczak et al .lll998l: 
Diamond & Kemball 1999; Bains et al. 2003; Vlemmin gs et alj 
l2005h . In addition, although very successful in simulating mass 
loss in carbon stars, frequency-dependent hydrodynamic mod- 
els are not able to reproduce the obs erved mass-los s rates in M- 
(C/0<1) and S-type (C/0~1) stars ('Woitkel l2006l) . unless spe- 
cial conditions are assumed (Hofner 2008). The observations of 
clumps and deviations from spherical symmetry, together with 
the difficulties to reproduce the observed mass-loss rates in M- 
and S-type stars, indicate that our current picture of mass loss on 
the AGB is perhaps too simple or lacking crucial ingredients. 

Late in the evolution on the AGB, the transition from 
(in most cases) a spherica lly symmetric CSE to an asymmet- 
ric pl anetary nebula (PN) (IZuckerma n & AUerlll986l Ide Marcol 
I2OO9I). where bipolar and elliptical mo rphologies are common 
(iMeixner et al.lll999l: lUeta et al.ll2000l) . continues to be a puz- 
zle. Several ideas o n how and when these features eme rge exist: 
binary interaction (iMorris IIT9871 iHuggins et al.ll2009l). interac- 
tion with a planet or a brown dwarf (e.g.,'Nordhaus & Blackman] 
2Q0©, or magnetic fields (e.g., Garcia-Segura et al. 2005). So 
far it is not clear which model (or combination of models) gives 
the more satisfactory explanation. Observations of remarkably 
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spherically symmetric detached shells around a handful of car- 
bon stars add to this mystery to some degree (e.g., Olofsson et al. 
1996, Gonzalez-Delgado et al. 2001, 2003, Maerckerel al. 2010, 
Olofsson et al. 2010). The general consensus is that the detached 
shells are formed as a consequence of a substantially increased 
mass-loss rate during a brief period, possibly following a ther- 
mal pulse. However, the details of the formation process, why 
they are only seen around carbon stars, or how their formation 
correlates with the potential onset of asymmetries toward the end 
of the AGB, is not clear 

Imaging the circumstellar envelope (CSE, gas or dust com- 
ponent) will give us general information on the symmetry of the 
mass loss, while comparing images at different wavelengths tells 
us something about the interaction between the two components. 
Both the general symmetry and the interaction are important for 
our understanding of how mass is expelled from the star and 
of the interaction further out in the CSE. Detailed images of 
AGB CSEs are so far available only for a limited number of (in 
most cases) high-mass-loss-rate sources. This makes it impossi- 
ble to know whether the asymmetries sometimes observed are 
a general feature of the AGB evolution or only present in some 
objects and are created only in combination with, for instance, 
a binary companion or a strong magnetic field. Problems with 
imaging the CSE are (for the gas component) insufficient spatial 
resolution at radio wavelengths, and (for the dust component) at 
shorter wavelengths, the very large star-to-CSE brightness ratio. 
One possible solution for the latter problem is provided by imag- 
ing polarimetry at optical and IR wavelengths. This is a compar- 
atively simple technique that requires relatively little telescope 
time, making it ideal to study large samples to be able to draw 
more general conclusions. 

PolCor (Polarimeter and Coronograph) is a new combined 
imager, polarimeter, and coronagraph that provides sharp images 
(resolution down to 0'.'2), and a well-defined point-spread func- 
tion (PSF) resulting in a high image contrast. In this paper we 
present a preliminary study to investigate its capability to image 
CSEs around AGB stars. The puipose is to study their structure 
and dynamical evolution. We have observed three AGB stars us- 
ing the PolCor instrument: the S-type binary AGB star W Aql, 
and the two detached shell sources DR Ser and U Cam. In Sect.|2] 
the imaging technique used in this work and previous investiga- 
tions are briefly discussed. In Sect. |3] the observations and the 
data reduction are described. In Sect. |4]the observed sources are 
presented. In Sect.|5]the analysis is outlined. In Sect.|6]the results 
are presented and later discussed in Sect. |7] Finally, a summary 
is given and conclusions are drawn in Sect. |8] The PolCor instru- 
ment and its performance are presented in Appendix lAl 

2. Imaging in polarized light 

When light is scattered by dust particles it becomes polarized. 
The intensity and polarization of the scattered light can be used 
to determine the properties of the dust and to map the dust distri- 
bution. The polarization degree is highest when the direction of 
the incident radiation is perpendicular to the viewing angle, but it 
will also depend on the wavelength, the grain size, and the grain 
composition. Zubko & Laor (200^ performed detailed calcula- 
tions of the spectral polarization properties of optically thin and 
thick dust in different geometries. They found that the degree of 
linear polarization in the optical and near-infrared is high (a;80% 
at 90° scattering angle) for wavelengths shorter than 0. 1 fim, de- 
creases up to about 0.2 yum, stays constant at around or below 
40% and starts increasing again long-ward of 1 //m. Only light 
that is scattered at an angle of 90° will be effectively polarized 



and polarized scattered stellar light hence probes the distribution 
of the dust in the plane of the sky. 

Imaging polarimetry has previously been proven to be a 
suitable technique for mapping structures in the close circum - 
stellar environment around post- AGB stars (iGledhilll l2005bh . 
iGledhill et alJ Jlool did ground-based near-infrared imaging 
polarimetry of 16 protoplanetary nebulae (PPNe). They found 
that a large majority of their objects were extended in the 
polarized-intensity images (as compared to the total-intensity 
images) showing that the objects are surroun ded by dusty en- 
velopes. This work was then followed up by IGledhilll (l2005al) 
where 24 additional sources were observed. Also here polariza- 
tion was detected in most of the observed sources. They sug- 
gest that the sources can be divided into objects with an opti- 
cally thick disk (probably due to binary interaction) resulting 
in a bip olar morphology, or optically thin dust shells. IGledhilll 
(I2005ah found maximum polarization degrees up to 60-70% in 
the bipolar objects. High spatial resolution observations using 
NICMOS on the HST have also been used to investigate the mor- 
phologies of dust envelopes around PPNe (lUeta et al.l2005h . and 
Su et al. (2003) found even higher polarization degrees around 
bipolar post- AGB stars. 

For stars on the AGB imaging polarimetry has been used 
when analyzing the large detached shells around ca rbon stars 
dGonzalez Delgado et all 120031: iMaercker et al.ll2010h . To esti- 
mate the widths and radii of the detached shells as accurately as 
possible is important for a better insight into how the shells are 
formed. It is also important in order to investigate whether the 
suggested He-shell flash scenario for the formation can be con- 
firmed, and for getting a better understanding of the mass-loss 
processes during the thermally-pulsing AGB phase (TP-AGB). 
Observations of the polarized light from the circu mstellar en- 
vironment of R Sc l ('Gonzalez Delga do et al.l 2003 ) and U Ant 
dGonzalez Delgado et al. 2003; Maerc ker et al.l 2010l) have pro- 
vided the possibility to measure the diameter and width of the 
dust shells with unprecedented accuracy. 

The advantage of observing scattered stellar Ught in the opti- 
cal in order to map circumstellar structure, especially compared 
to single-dish radio observations, is the high spatial resolution 
along with high-sensitivity detectors. In the optical, the star-to- 
CSE flux ratio is ~ lO'*, requiring a very large dynamical range of 
the detector in or der to image the weak emission from the CSE 
( lUeta et a"ni2000h . Assuming that the direct light from the cen- 
tral star is essentially unpolarized, it should disappear in images 
of polarized light, reducing the required dynamical range of the 
detector However, already a small amount of polarization will 
leave a significant signature of the stellar PSF also in the polar- 
ized images, and it is therefore sometimes necessary to also use 
a coronograph. 



3. Observations and data reduction 

3.1. Observations 

The observations were carried out during three different obser- 
vation runs at the 2.5 m Nordic Optical Telescope (NOT) on La 
Palma, The Canary Islands, using the PolCor instrument. In June 
2006 we used a prototype of PolCor based on the same piinci- 
ples as the final instrument (described in detail in Appendix lAli. 
The observations are summarized in Table [1] 

During each observation run we observed a number of po- 
larizat ion standard stars c hosen from the lists of Turnsh ek et al.l 
( 119901) and lHeilesI ( |2000|) . The results of the calibration for the 
observations in July 2008 are shown in Fig. [1] The left panel 
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Table 2. Variable type, spectral type, period, P, distance, D, 
luminosity, L*, and present mass-loss rate, M, of the observed 



Pd_stnd [%] 



100 150 
V'-Stnd [degrees] 



Fig. 1. Calibration observations with PolCor in the V band for 
the July 2008 observation run. The left panel shows the degree 
of polarization of the calibration stars at the time of observa- 
tion compared to the literature values. The right panel shows the 
relation between the observed polarization angle and that listed 
for the calibration stars. The off-set angle (22°8) is due to the 
mechanical attachment of the instrument to the telescope. 

Table 1. The PolCor/NOT observations. Dates of the observa- 
tions, the filters and masks used, the exposure time, t, and the 
effective seeing in the final total-intensity images, are given. 



Source 


Date 


Filter 


Mask 


t 


Seeing 








["] 


[h] 


["] 


W Aql 


2006-06-10 


V 


3 


1 


0.8 


W Aql 


2006-06-10 


V 


open 


1 


0.7 


W Aql 


2008-07-05 


R 


6 


0.5 


0.7 


W Aql 


2008-07-05 


R 


open 


1 


0.7 


UCam 


2007-10-06 


V 


3 


1 


0.8 


DRSer 


2006-06-10 


V 


3 


1 


0.7 



shows the degree of polarization, and the right displays the po- 
larization angle. The formal error bars of the measurements are 
smaller than the symbols. The scatter in the left panel is most 
likely due to time variations (most of the literature values are 
old) and/or it reflects differences in the effective wavelength. In 
the right panel a close relation between the observed polarization 
angle and that listed for the calibration stars is shown. The off-set 
angle (22°8) is due to the mechanical attachment of the instru- 
ment to the telescope. The polarization of the NOT is negligible 
(A. Djupvik, NOT Senior Staff Astronomer, private communica- 
tion). 



3.2. Data reduction 

The 'Lucky imaging' technique was used to improve the image 
quality. This means that only the sharpest frames were selected 
and used to produce the final composite image. The frame rate 
of the observations was 10s"'. Both the sharpness and the image 
motion were calculated for each frame. The sharpness is defined 
as the percentage of the total light inside a O'.'55x0'.'55-box cen- 
tered on the brightness peak. The image motion is determined by 
the location of the brightness peak compared to an average value. 
There is a trade-off between sharpness and depth and an appro- 
priate sharpness acceptance level was chosen for each observed 
source and wavelength band. Shift-and-add procedure was per- 
formed only for the accepted frames. At a seeing of 0'.'? just 
shift-and-add typically improves the measured seeing to Of! 5. If 
another 50% of the frames are excluded, the seeing is further im- 
proved to about 0'.'4. The seeing at the different observation runs 
was slightly worse than 0'.'7. The effective seeing (Full Width at 



sources. 


Source 


Var. 


Spec. 


P 


D 




M ' 




type 


type 


[days] 


[pc] 


[Lo] 


[Moyr-'] 


W Aql 


M 


S 


490 


230- 


6800 


22x10-' 


HD204628^^ 




FO 










UCam 


SRb 


C 


400 


430-' 


7000 


2.0x10-^ 


DR Ser 


Lb 


C 




760^ 


4000 


0.3x10-^ 



' From iRamstedt et al.l ( 120091) (W Aql) and ISchoier et all ( 120051) 

(DR Ser and U Cam). 

- From period-luminosity relation (IWhitelock et al.|[l994t) . 

^ Standard star used in the analysis of W Aql. 

■* From period-luminosity relation (iKnapp et al.ll2003]) . 
' From adopting = 4000 Lq. 



Half Maximum, FWHM) of the final images after shift-and-add 
and sharpness-selection procedures is given in Table [1] Further 
details of the performance of the PolCor instrument and the data 
reduction are given in Appendix lAl 

After frame-selection and co-adding procedures we get im- 
ages at the four polarizer positions plus a dark image. All im- 
ages are quasi-simultaneous. Changing atmospheric conditions 
as well as instrumental drifts cancel when defining the Stokes 
parameters, I, Q and U: 



/ (A0-I-A45-HA90-HA135) 

A X I U 



Q = A0-A90 
U = A45 - AOS, 



(1) 

(2) 
(3) 



where the image at a polarizer position of degrees is denoted as 
AO, etc., and the dark image is denoted as D. The sky emission 
is defined as the median value of corner regions away from the 
PSF of the star and subtracted for I, Q and U. Then we calculate: 



p = VG- + u~ 

P 

(U 

il/ — 0.5 atan — 

\Q 



(4) 
(5) 

(6) 



where P is the polarized intensity, P^ the polarization degree and 
i/r the polarization angle (counted from north toward east). 



4. Observed sources 

4.1. The S-type AGB star W Aql 

W Aql is an S-type Mira variable with a period of 490 days 
dKukarkin et aLlfigTlh . This correspond s to a distance of 230 p c 
using the period-luminosity relation of IWhitelock et al.l (Il994l) . 
No Hipparcos parallax has been measured. It w as first identi- 
fied as a spectroscopic binary by iHerbid (Il965b . He suggested 
that the companion is an F5 or F8 star and that the separa- 
tion is less than 0'.'8. Figure |2] shows the high-resolution im- 
age of W Aql observed with the Advanced Camera for Surveys 
(ACS) on the Hubble Space Telescope (HST) in High Resolution 
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2 1 0-1-2 

east offset [arcsec] 

Fig. 2. The high-resolution HST B-band image of W Aql. The 
binary system is clearly resolved. North is up and east is left. 

Channel (HRC) mode using the F435W filteiQ. The spatial reso- 
lution of the image is (X' 12. The image was downloaded from the 
HST archive^ The binary pair is clearly resolved and the stars 
are separated by 0'.'46. This corresponds to x 110 AU assuming 
a distance of 230 pc. The inclination of the orbit is unknown and 
the orientation of the sources relative to each other is therefore 
also unknown. In Fig. |2] they appear to be oriented in a north 
east-south west direction, but this could very well be a projection 
effect. 1 10 AU is therefore the minimum separation. This corre- 
sponds to a minimum orbital peiiod of x 630 yrs if a total mass 
of the binary system of 3 Mq is assumed. The mass-loss rate of 
W Aql is estimated to be 2.2 x 10'^ Mp yr~' fro m observations 
of CO radio line emission (iRamstedt et al.ll2009h . 

The emission from the star and its surrou ndings has been 
exten sively observed, both in th e infrared (e .g., [Tevousian et al.l 
l2004t l^nkata Raman & Anan darad l2008h and i n the radio 

ISchoier et a 



regime (e.g.. IDecin et al.ll2008t ISchoier et al.ll2009h . This is the 
first time measureme nts of the po l arized light around W Aql 
have been published. iTatebe et al.l (l2006l) observed W Aql at 
ll.lSyum with the three-element UC Berkely Infrared Spatial 
Interferometer (ISI) and produced one-dimensional profiles of 
the dust emission within the nearest O'.'S of the star. They found 
that the dust emission was enhanced on the east side of the star 
relative to the west side. They also found that the intensity de- 
cline on the east side was more gradual. They interpreted their 
results as indicative of a dust shell that has been expelled ap- 
proximately 35 yrs prior to their observations. 

4.2. The detached shell sources U Cam and DR Ser 

In addition to the 'normal' CSEs observed due to the continuous 
mass loss from AGB stars, large and geometrically thin shells are 
found to exist around a small number of carbon stars. These de- 
tached shells were first detected by Olofsson et al. (1988), who 
found wide double-peaked profiles in the CO line emission from 
the carbon stars U Ant and S Set. They have been suggested 
to be a consequence of the increase in luminosity during the 



' HST Proposal 10185 PI: Raghvendra Sahai 
Title: When does Bipolarity Impose itself on the Extreme Mass 
Outflows from AGB stars? An ACS SNAPshot Survey 

- Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 



He-shell flash/thermal pulse (lOlofsson et al.|[T990l) . A more de- 
tailed study of the properties of the detached shells was done by 
Ipiofsson et al. ( 1996) using maps of the CO line emission. For 
the stars R Scl, U Ant, S Set, V644 Sco and TT Cyg they found 
spherically symmetric shells with radii between 1 - 5 x 10'^ cm, 
expansion velocities of 13-20 kms"', and ages of 1 - 10 x lO-' 
years. The shells are thin {AR/R < 0. 1 ) and high-resolution maps 
observed with the IRAM Plateau de Bure Interferometer (PdBI) 
confirm their remarkable spherical symmetry (Lindqvist et al. 
1999; Olofsson et al. 2000). The detached shells around R Scl 
and U Ant were observed in scattered light in direct imaging 
mode and in polarization, showing the distribution of the dust 
and the gas in the shells (Gonzalez-Delgado et al. 2001, 2003; 
Maercker et al. 2010). The detached shells around R Scl and 
U Cam were observed in dust scattered light with the HST, show- 
ing a significant amount of dumpiness in the shells (Olofsson et 
al. 2010). 

The stellar parameters for U Cam and DR Ser are given in 
Table m These shell s are relatively small (also in spatial scale). 
ISchoier et al.l (l2005h performed detailed radiative transfer mod- 
eling of both the molecular linr emission and the spectral en- 
ergy distribution (SED) of the seven currently known detached- 
shell sources. The present-day stellar mass-loss rates and the 
masses of the detached shells were estimated. For U Cam the 
shell mass and mass-loss rate are estimated to be 1 x 10""' Mq 
and 2 x 10"^ Mq yr ', respectively. For DR Ser the correspond- 
ing values are found to be 1 x 10"^ M© and 3x lO^^^Moyr"', 
respectively. The expansion velocities of the shells are deter- 
mined through the molecular line emission and are found to be 
23kms ' and 20kms"' for U Cam and DR Ser, respectively. 
The shell around U Cam has a radius of 7?sh =4.7 x lO'^cm 
(measured from interferometric CO observations), while the 
shell around DR Ser has R^u = 8 x 10"'cm (found from fitting 
on-source CO radio line spectra). 



5. Analysis 

5. 1 . Circumstellar structure 

The images of polarized light can be used to determine the struc- 
ture and the physical extent of the circumstellar dust envelope 
and of the detached dust shells (i.e. the shell radii and widths). 
The PolCor coronographic mask reduces the brightness of the 
central stars by a factor of 100, making it possible to detect very 
faint circumstellar light. 

For both the polarized- and total-intensity images, the stellar 
PSF was approximated by fitting a polynomial (4th degree) to an 
azimuthally averaged radial profile (AARP) in loglog-scale. The 
PSF was then subtracted from the original image to emphasize 
weak features. New AARPs were then calculated covering dif- 
ferent sections of the images to find the extent and structure of 
the different circumstellar features. Distortions due to the tele- 
scope spiders were avoided by averaging over appropriate an- 
gles. 

To determine /?sh and A/?sh of the detached shells, the PSF- 
subtracted AARPs of the polarized intensity were fitted, assum- 
ing isotropic scattering by dust grains. The calculation also as- 
sumes that the dust follows a Gaussian radial density distribu- 
tion in a shell of radius R^t, and with a FWHM of A/?sh (see 
iMaercker et ani2010l:loiofsson et al.ll2010i for details). 
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5.2. Calculating the dust mass 

Following the procedure of iGonzalez Delgado et al.l (1200 ll) . we 
have performed a simple analysis to calculate the total dust mass, 
Md, of the observed structures. If an optically thin dust envelope 
is assumed, the scattered flux, F^^, is given by integrating the 
product of the stellar flux (at the observed distance from the star), 
Fi^lAnR^, the number of scatterers, A^sc, and the scattering cross 
section, cTsc, over the observed wavelength range: 

^'sc = J fr,,F^dA. (7) 

and were measured in the total-intensity images. 
was estimated by fitting a Gaussian to the star and summing all 
counts within 3cr. This was then corrected for the 5-magnitude 
dampening of the mask. The scattering cross section for spheri- 
cal grains is given by 

O-sc = GscTTfl^, (8) 

where Qsc is the scattering efliciency and a is the grain radius, 
for simplicity assumed to be constant at 0. 1 jjm. 

The ratio between Fsc and gives A^sc through Eq. |7] The 
mass of each grain, mj, is given by the volume and density of 
the dust grains. For the silicate and carbon grains we assumed 
typical values for the grain density of 3 and 2 g cm"^, respec- 
tively (e.g.. lSui3[T999[l2000[ for silicate and carbon grains, re- 
spectively). The total dust mass, Md, is then given by A^sc x 'Wd- 
For W Aql, a stronomical silicates with optical constants from 
iDraind ([T98l were used to calculate Qsc- For the c arbon stars , 
U Cam and DR Ser, amorphous carbon grains from ISuh] (I2OOOI) 
were used. The estimated total dust mass depends strongly on the 
assumed grain size. Mi oc aj gsc, where Qsc cfi, and 2 < /? < 3 
in the wavelength range of our observations. The assumption of 
a constant grain size further adds to the uncertainty in the mass 
estimate. In addition, there are a number of uncertainties and 
simplifications in this method and the estimates are therefore 
order-of-magnitude estimates. 

6. Results 

6. 1 . The dusty environment around W Aql 
6.1.1. Circumstellar structure 

To study the circumstellar dust distribution around W Aql the 
R-band images acquired with the 6" coronographic mask were 
used. Already in the total-intensity image (Fig.|3] left) the bright- 
ness distribution around W Aql appears asymmetric. The scat- 
tered light is more intense to the south-west (SW) where the 
emission extends out to about 10" from the center, compared 
to about 5" in the north and the east directions. 

The polarization angle and the degree of polarization are 
shown as polarization vectors in Fig. ^(righf). The vectors are 
overlaid on an image showing the polarized intensity. Figure [3] 
(right) shows that the polarizing dust is distributed all around the 
star, however, the SW enhancement appears clearly. The image 
has been smoothed by a 3x3-pixel Gaussian to reduce the noise. 
The integrated polarization degree across the image (correspond- 
ing to what would be measured if the source was unresolved) is 
about 10%. In the SW quadrant the mean polarization degree is 
around 20% across the feature. The maximum polarization de- 
gree is found to be just above 40% in the SW part of the image. 

Figure|4](/ar left) shows the polarized intensity within 10" of 
W Aql. To investigate the distribution and extent of the SW 



feature, the image was divided into four quadrants: north-east 
(NE), north-west (NW), south-east (SE), and south-west (SW). 
By calculating AARPs of the different quadrants, we can com- 
pare the brightness distribution across the image (Fig. |4] middle 
right). The image is clearly brighter in the SW, while the other 
three quadrants look similar. A 4th degree polynomial (dashed 
line. Fig. |4l middle right) is well-fitted to the log-log AARP of 
the three weaker-emission quadrants (solid line. Fig. |4] middle 
right), and by subtracting the fit from the polarized-intensity im- 
age, the location and extent of the SW asymmetry was found 
(Fig.m/ar right). The area close to the mask (<4'.'5 offset) and 
the outer parts (>16"offset) were not taken into account in the 
fitting. The SW brightness enhancement seems to start already 
at the edge of the mask and it is nearly constant out to about 7". 
It then declines until it disappears at approximately 12". Figure|4] 
(middle left) shows the image after the fit to the weaker-emission 
quadrants has been subtracted. 

W Aql was also observed without the coronograph both in 
the R- and V-band. The polarized-intensity images are shown in 
Fig- El overlaid by polarization vectors. In the R-band image the 
star was saturated and any structures in the very inner parts of 
the image can therefore not be trusted. This also explains the 
read-out tracks seen in Fig. ^(left). The polarized intensity of 
the V-band image is shown in Fig. |5] (right). To interpret the 
very inner structure is not straightforward as the secondary star 
contributes significantly in the V-band and the binary system is 
not resolved. Both images in Fig. |5] show the same polarization 
pattern as can be seen in Figs|3] and|4]and the SW feature is also 
clearly confirmed. 

6.1 .2. The dust mass of the SW feature 

The amount of dust in the SW brightness enhancement was 
calculated assuming optically thin dust scattering (Sect. 15.2b . 
By measuring the ratio between the scattered flux and the stel- 
lar flux in the total-intensity image, N^c was obtained accord- 
ing to Eqn. I2] The scattered flux was estimated by adding all 
counts in the SW quadrant of the image from 4" (to avoid dis- 
tortions due to the edge of the mask) out to 12". The stellar 
flux, Fi,, was estimated as described in Sect. 15.21 The flux ra- 
tio was found to be Ff^JF-^ = 3 x lO"-', and the dust mass was 
found to b e Md ~ 1 x 10'^ M(T). Assuming a dust-to-gas ratio of 
1 . 1 X 10"^ dRamstedt et alJ20"09l) . this corresponds to a total mass 
of 10-^ Mq. 

6.1 .3. A close-up view of W Aql 

The coronographic masks used during the observations are not 
entirely opaque, but only dampen the light by 5 magnitudes. This 
opens up the opportunity to investigate the close circumstellar 
environment as seen through the mask (without saturation prob- 
lems as in Fig. |5] left). W Aql was observed with the 6" mask 
in the R-band and with the 3" mask in the V-band. The V-band 
image attained with the smaller mask is distorted by spill-over 
from the sides and by emission from the secondary star and 
therefore not shown. Figure |6] shows the contour maps of the 
total (left) and polarized (right) R-band intensity seen through 
the 6" mask. The total-intensity image clearly shows the pri- 
mary star and that its shape is not distorted by the mask. The 
secondary star is not visible in the R-band, instead the image 
is dominated by the emission coming from the AGB star. In the 
polarized-intensity image the close environment seems stretched 
in the NE-SW direction toward the SW asymmetry seen in the 
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Fig. 3. Left: The total-intensity image of W Aql. Right: The polarized intensity around W Aql overlaid by polarization vectors 
showing the polarization degree and angle. North is up and east is left. 




east offset [arcsec] east offset [arcsec] log(offset) offset [arcsec] 



Fig. 4. The R-band polarized-intensity image of the circumstellar envelope around W Aql. North is up and east is left. Far left: The 
polarized intensity, P, of the inner 10". Middle left: The polarized-intensity image after the fit to the three weaker-emission quadrants 
(SE, NE, N W, see text for further explanation) has been subtracted. Middle right: The log-log AARPs of the polarized intensity over 
the different quadrants. The y-axis is normalized to the polarized intensity "of the star", P», i.e., the polarized intensity within the 
3cr area of a Gaussian fitted to the star (see Sect. l5.2l for details on how the size of the star was estimated). The vertical dashed lines 
indicate the most reliable region of the image, and this is also the area used for the fit. The areas close to the coronographic disk, and 
close to the edge of the image, are excluded since they are less reliable. The thin dotted line shows the AARP over the NW quadrant. 
The thin dashed line shows the AARP over the SE quadrant. The thin solid line shows the AARP over the NE quadrant. These three 
quadrants look very similar in the reliable region and the thick solid line shows the average of all three. The thick dash-dotted line 
shows the SW quadrant, which is clearly brighter in this area of the image. The thick dashed line shows the polynomial fit to the 
average over the three weaker-emission quadrants. Far right: The AARP over the SW quadrant after the polynomial fit to the three 
weaker-emission quadrants has been subtracted. 



wide-field images, and there is an indication of a bipolar struc- 
ture. The detached shell sources were also inspected through the 
mask, and although U Cam was slightly aspherical, neither has a 
similar structure to that found in Fig. |6] Further imaging obser- 
vations of W Aql would be necessary in order to firmly confirm 
this tentative structure. 



6.2. The detached shells around U Cam and DR Ser 
6.2.1. Circumstellar structure 

In images of polarized light, detached shells appear as ring-like 
structures with no intensity at small (projected) distances from 
the star (e.g., Fig.|7] upper row, far left). In the total-intensity im- 
ages, however, the scattered light may significantly contribute to 
the intensity measured at small radii, depending on the forward 
scattering efficiency of the grains. This leads to images that have 
a more disk-like appearance, possibly with some degree of 'limb 
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Fig. 5. Left: The polarized intensity of W Aql in the R-band observed without the coronographic mask. The star is saturated in 
this image and any structure seen in the very inner parts can therefore not be trusted. This also explains the read-out tracks in 
the image. Right: The polarized intensity of W Aql in the V-band observed without the coronographic mask. The secondary star 
contributes significantly in the V-band, but the binary pair is not resolved. Both images are overlaid by polarization vectors showing 
the polarization degree and angle. North is up and east is left. 




eastoffset[arcsec] eastoffset[arcsec] 

Fig. 6. Left: Contour map of the R-band total intensity as seen through the 6"mask. The contours are drawn in log-scale at 10%-steps 
relative to the maximum intensity. The outermost contour of the central structure corresponds to an 80% decrease in the intensity. 
Right: Contour map of the R-band polarized intensity as seen through the 6"mask. The contours are drawn at 10%-steps relative 
to the maximum intensity. The outermost contour of the central structure corresponds to a 70% decrease in the intensity. In both 
images, the cross marks the location of the star and the thick line marks the edge of the mask. North is up and east is left. 



brightening' (lOlofsson et al.l l20To[ see also e.g.. Fig. [T] lower 
row, far left). 

Figures |7] and [8] show the images of the polarized intensity, 
P {upper row, left), and the corresponding AARPs {upper row, 
right) of U Cam and DR Ser, respectively. The polarized inten- 
sity in the shells, Psh, and the corresponding AARPs are given 
in the/flr left and middle right panels. Masks have been placed 
over the inner parts of the images where the data is not reliable in 
the middle left panels. When creating the AARPs the telescope 
spider was avoided by averaging only over the position angles 
between 73° - 115° and 250° - 315° (U Cam), and 15° - 90° 
and 210° - 260° (DR Ser) and the images were smoothed us- 
ing a Gaussian kernel with cr=0'.'12 and 0'.'092 for U Cam and 
DR Ser, respectively. The polarized intensity in the shells {mid- 
dle left) was attained by subtracting a fit to the PSF in the middle 
right panels (dashed line), as described in Sect. 15.11 When fit- 
ting the PSF, the inner regions that are significantly affected by 
the mask (radii < 2"), and the regions where the shells are lo- 



cated (indicated by the vertical dashed lines in Figs [T] and [8]l, 
were avoided. The radii and widths of the shells were estimated 
following the procedure of Maercker et al. (2010) and Olofsson 
et al. (2010) as described in Sect. 15.1 1 The results are given in 
Table H 



6.2.2. The dust masses of the detached shells 

The images of the total intensity / towards U Cam and DR Ser 
are shown in Figs |7] and |8] {lower row, left) together with the 
respective AARPs {lower row, right). The AARPs were created 
averaging over the same ranges of position angles as for the po- 
larized images. The same regions were excluded when fitting the 
stellar PSF as for the polarized images (indicated by the vertical 
dashed lines). Masks cover the inner regions of the images where 
the data is not reliable in the upper right panels. To determine 
the circumstellar flux, F^^, a constant intensity distribution out 
to radius /?sh was assumed. At /?sh, the assumed intensity profile 
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Table 3. Results of the fits to the polarized and total intensity 
AARPs for U Cam and DR Ser 



7.1.2. Possible shaping agents 



Source 


["] 


AR.k ["] 




Mi [Mo] 


UCam 


7.9 


0.9 


5 X 10-^ 


5 X 10-' 


DRSer 


7.6 


1.2 


9 X 10-* 


2 X 10-^ 



then follows a Gaussian decline with a FWHM of AR^h- Rah and 
A/?sh were determined from the polarized images as described 
above. The assumed radial profiles are shown as dashed lines in 
the lower middle right panels of Figs [T] and |8] The stellar fluxes 
were estimated as described in Sect. 15.21 The derived Fsc/^^* 
ratios and dust masses are given in Table |3] 

Using previous results from radiative modeling of the molec- 
ular emission (see Lindqvist et al. 1999; Schoier et al. 2005, for 
U Cam and DR Ser, respectively) the dust-to-gas ratio dur- 
ing the creation of the shells can be estimated. For the shells 
around U Cam and DR Ser the gas masses are estimated to be 
1 X IO^^^Mq. This leads to a dust-to-gas ratio of 5 x lO""* for 
U Cam and 2 x 10"^ for DR Ser 



7. Discussion 

7.1. The dusty environment around W Aql 

7.1 .1 . Tine circumstellar dust distribution 

The estimated dust mass together with the polarized-intensity 
images show that there is a significant amount of dust around 
W Aql (Sect.l4m). The dust is found ah around the star (Figs|3]|5] 
and |6|, but the images (Figs [3] H) |5] and ^ also clearly show 
that the dust is distributed asymmetrically. The brightness of the 
scattered light is enhanced on the south-west side of the star. The 
degree of polarization of the light from the SW feature is higher 
than on the other sides of the star The polarization degree in the 
SW is in Une with what can be expected when light at 0.6 jum 
is scattered at a 90° angle by optically thin dust, but higher than 
what would be expected for optically thick dust (IZubko & Laorl 
120001) . The increased polarization degree further supports the 
suggestion that there is an increase in the amount of dust on the 
SW side of the star 

An asymmetr ic dust di stribution around W Aql was dis- 
cussed already bv iTatebe et a l. (2006). They find that the asym- 
metry probably has been stable for a "moderately long time 
span", and at least for 35 yrs. From this they conclude that a 
companion is an unlikely source of the asymmetry, since 35 yrs 
is of the order of the orbital period required for a close com- 
panion to be responsible for the shaping. The asymmetry in the 
dust distribution is apparent even far from the star with consis- 
tently more dust on the east side. They conclude that this may be 
due to convection and turbulence, or magnetic eff'ects, causing 
material to be ejected preferentially on one side of the star In 
this context i t is im portant to point out that the observations of 
iTatebe et al.l (l2006b only probe out to 500 milliarcseconds and a 
direct comparison is therefore difficult. However, from the ob- 
servations presented in this work it is clear that the asymmetry 
reaches further out. 

Some asymmetry or structure can also be distinguished in 
the CO ra dio line profiles, whic h trace the circumstellar gas dis- 
tribution dRamstedt et al.l 120091) . The CO lines appears slightly 
brighter on the blue-shifted side, i.e., from the front side of the 
CSE where gas is moving towards us. This could possibly indi- 
cate an increase also in the gas density. 



According to the classification of binary interaction by ISokej 
(Il997h . W Aql probably falls within the category of wide bina- 
ries (100 AU < fl < 1000 AU, where a is the orbital separation). 
The most likely outcome of a wide binary when evolving into 
a PN is that the nebula is asymmetric. The SW asymmetry in 
the CSE of W Aql appears to be aligned with the binary orbit 
(Fig.|2]i and the binary separation is at least 1 10 AU. The orbit is 
not known and it is possible that the seen alignment is merely a 
projection effect with one star placed in front of the other The bi- 
nary separation of 1 10 AU is thus a lower limit. This uncertainty 
and the lack of theoretical models with large binary separation 
(> 70 AU), makes it is unclear whether binary interaction could 
cause the observed asymmetry. 

The well-studied archetypic a l mira star o Ceti (see, e.g., 



work bv iKarovska et al.1 1 19971 120051: iMatthews & Karovskal 



2006: M atthews etalJ l2008r and references therein) serves as a 
good comparison as it is also a binary star and the binary sepa- 
ration (fl ~ 100 AU) is of the same order as, or smaller than that 
of W Aql. The Mira AB system ha s been observed on different 
scales and at different wavelengths. [Ireland et al.l (|2007l) studied 
the close circumstellar environment at mid-infrared wavelengths 
and found that Mira B is a main-sequence star (the data is found 
consistent with a 0.7 Mq K5 dwarf) surrounded by a 10 AU ac- 
cretion disk. The nature of Mira B is however still a matter of 
debate an d there are strong indication s that it is a white dwarf 
(see e.g.. ISokoloski & Bildstenll2010l and references therein). 
Mira B is found to have an accretion radius of 3 2 Mb /M(T)AU 
(assuming Bondi-Hoyle-Lyttleton accretion, e.g., lEdgaill2004h . 
If we perform the same calculation for W Aql (see Sect. 3.4 in 
the paper by jlreland et al. 2Q(}2), although noting that the param- 
eters of this system are even less known than for Mira AB, the 
accretion radius of W Aql B is less than 6 Mb /Mq AU, making 
accretion from W Aql A to W Aql B less likely. 

Mira A is kno wn to have a bipola r outflow (studied in 
molecular emission; Planesas et al .Il990allbl:ljo"sselin et al.l20"00l 
Fong et al. 2006) on approximately the same scale as the asym- 
metry in W Aql. Also in a well-studied case like Mira, under- 
standi ng the origin of circu mstellar structure is not straightfor- 
ward. iJosselin et al.l (|2000|) suggest that asymmetric mass-loss 
due to e.g. non-radial pulsations, giant convection cells, and/or 
magnetic spots on the surface provide the most likely expla- 
nati on for the asymm etries seen in the CO emission from the 
star lFong et al.l(l2006h find that their observations support a sce- 
nario where the structures are caused by interaction between the 
molecular gas and a rotating disk. 

Observations required to investigate the kinematical proper- 
ties of the asymmetry seen in W Aql are not available at this 
point. Interferometric spectroscopy should be performed in or- 
der to determine the velocity field and to investigate whether 
the structures seen can be an earlier version of the outflow seen 
from Mira A. The CO emission from W Aql has been ob- 
served at the PdB interferom eter as part of the COSAS program 
(ICastro-Carrizo et a"Lll2010b . The results are planned for a future 
publication. 

In order to properly evaluate the origin of the asymmetrical 
dust distribution around W Aql, a hydrodynamical simulation 
using the known parameters of this system needs to be done (for 
models of a separation up to 70 AU see de Val::Borro et al. 20091). 
The possible shaping due to a magnetic field cannot be evaluated 
at this point as it has not yet been measured around this star. 
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Fig. 7. The detached shell around U Cam. North is up and east is left. The upper row shows the polarized intensity images. The 
lower row shows the total-intensity images. Upper row, far left: The polarized intensity P. Upper row, middle left: The polarized 
intensity in the shell Psh- Upper row, middle right: The AARP of the polarized intensity normalized to the stellar intensity (same as 
in Fig. m. The dashed line shows the PSF-fit. The vertical dashed lines indicate regions not included in the fit, i.e., the inner part 
close to the mask and the region where the shell is located. The vertical dotted line shows the determined radius of the shell (see 
text for details). Upper row, far right: The AARP of the subtracted image. The lower row shows the same images but for the total 
unpolarized intensity / (far left and middle right), and the total unpolarized intensity in the shell Ish (middle left and far right). 



7.2. The detached shells around U Cam and DR Ser 

7.2.1 . Thie phiysical parameters of the shells 

The determined radius and width of the detached shell around 
U Cam (Table [3]l are in excellent agreement with the results of 
lOlofsson et al.l (uOlOh . who derived a shell radius of 7'.'7 and a 
width of (y.'6 based on HST images. At the distance of U Cam, 
the radius and width presented here coiTespond to 5 . 1 x lO"' cm 
and 5.8 x 10'^ cm, respectively. Measurements of the size of 
the CO shell based on data from the PdB i nterferometer re- 
sult in a radius of 7'.'3 dLindqvist et alJ[r999h . consistent with 
our results. However, the smaller size determined for the de- 
tached CO shell may indicate that the dust and gas have sepa- 
rated, similar t o the situation for the detached shell source U Ant 
dMaercker et al . 2010). The resulting d ust mass is a factor o f 
^2 higher than the mass determined by lOlofsson et al.l (1201 Oh . 
well within the uncertainties. The derived dust-to-gas mass ra- 
tio for U Cam (5 x 10""*) is lower than what has previously been 
found for carb on stars (2.5 x lO""* for sources with P < 500 days, 
iGroenewegen et al. 1998i) . 

Although the detached shell around D R Ser has previousl y 
been observed in CO spectral line emission (ISchoier et al.l2005l) . 
this is the first time the detached shell around DR Ser has been 
imaged. Thus, it is the first time that the size and the width 
of the dust shell around DR Ser have been measured directly. 
At the distance of DR Ser, the radius and width correspond to 
8.6 X 10'^ cm and 1 .4 x lO'*" cm, respectively. However, the large 
distance to the object (760 pc) makes it likely that the shell width 
given here is only an upper limit. Models of the molecular CO 
line emission give a radius of 8 x lO'*" cm, based only on the line 



shape. The dust mass in the shell derived from the PolCor dat a 
is a factor >5 lower than the estimate by ISchoier et al.l (l2005l) . 
well within the uncertainties (which are particularly large for the 
Schoier et al. result). The estimated dust-to-gas ratio during the 
formation of the shell around DR Ser (2 x lO ^^) is in good agree- 
ment with previous results for carbon stars (iGroenewegen et aU 
19981). 

For both U Cam and DR Ser the assumed AARPs of the total 
intensity differ significantly from the observed profiles, in par- 
ticular at small radii. The observed radial profiles are, however, 
very sensitive to the PSF subtraction, contributing significantly 
to the uncertainty in the determined dust masses. To some de- 
gree potential dumpiness in the shells and varying shell radii 
and widths will further affect the observed radial profiles. 

7.2.2. A connection to thermal pulses? 



Several previous investigations suggest a conne ction between 



the creation of detached shells and thermal pulses (Schoier et al. 


2005 


jMattsson et al.i2007£iMaercker et al.l2010tiOlofsson et al. 


2010 


). The previously measured CO expansion velocities to- 



gether with the extents of the shells measured in this work 
give upper limits for their ages of rucam~ 700 yr and ToRSer^ 
1400yr. The widths of the shells imply formation times of 
Arucam~ 100 yr and AroRSei * 200 yr. Based on the masses con- 
tained in the shells and the formation times, a simple calcu- 
lation implies mass-loss rates of ^2 x lO^^Moyr"' and x!5 x 
10 ''Moyr^' during the formation of the shells for U Cam 
and DR Ser, respectively. The present day mass-loss rates are 
2.0 X 10"^ Mq yr"' and 3.0 x 10"** Mq yr"', respectively. This is 
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Fig. 8. The detached shell around DR Ser. North is up and east is left. The upper row shows the polarized intensity images. The 
lower row shows the total-intensity images. Upper row, far left: The polarized intensity P. Upper row, middle left: The polarized 
intensity in the shell Psh- Upper row, middle right: The AARP of the polarized intensity normalized to the stellar intensity (same as 
in Fig. m. The dashed line shows the PSF-fit. The veitical dashed lines indicate regions not included in the fit, i.e., the inner part 
close to the mask and the region where the shell is located. The vertical dotted line shows the determined radius of the shell (see 
text for details). Upper row, far right: The AARP of the subtracted image. The lower row shows the same images but for the total 
unpolarized intensity / (far left and middle right), and the total unpolarized intensity in the shell (middle left and far right). The 
peak at 15"in the lower right plot is due to a background star. 



in line with the change in mass-loss rate during the thermal pulse 
re sponsible for the crea tion of a detached shell in the models 
of iMattsson et alJ (l2007h . If the shells are created due to a two- 
wind interaction scenario, the exact details of the interaction will 
significantly complicate a straightforward interpretation, and it 
must be stressed that the above estimates are order-of-magnitude 
estimates. We conclude that our results are consistent with pre- 
vious investigations and with the thermal-pulse-formation sce- 
nario. 



8. Summary and Conclusions 

We have investigated how the new imaging polarimeter and 
coronograph PolCor can be used to study the circumstellar du st 
distribution around AGB stars. In this preliminary study, ob- 
servations of the circumstellar structure around the S-type star 
W Aql and the two detached-shell sources, DR Ser and U Cam, 
were performed. Here we summarize our results and draw the 
following conclusions: 

- The images of W Aql show that the circumstellar dust dis- 
tribution is asymmetric, both on large (~10") and on smaller 
(~1") scales. The wide-field images show what appears to 
be a dust-density enhancement on the south-west side of the 
star 

- The polarization degree is found to be consistent with what 
could be expected when the incident light is scattered 90° by 
optically thin dust. 



- The dust mass of the SW feature around W Aql is estimated 
(assuming optically thin dust scattering) to be w 1x10"^ M©. 

- The close circumstellar environment around W Aql, as seen 
through the coronographic mask, exhibits an elongated, pos- 
sibly bipolar structure around the AGB star 

- Further observations to determine the kinematics and ob- 
tain information about possible magnetic forces, as well as 
hydrodynamical modeling to investigate the interaction be- 
tween the binary pair, should be performed in order to inves- 
tigate the cause of the asymmetric dust distribution around 
W Aql. 

- The detached shells around U Cam and DR Ser can be 
clearly seen in the polarized images. This is the first time the 
detached shell around DR Ser has been imaged. The radii 
and widths of the shells are determined and found to be con- 
sistent with previous results from imaging of CO radio line 
emission and HST images of dust-scattered light (U Cam) 
and from CO radio line shape modeling (DR Ser). For 
U Cam the radius and width are 5 x lO"" cm and 6x10'^ cm, 
respectively. For DR Ser they are 9x10'^ cm and 1 x 10'^ cm, 
respectively. 

- The total dust masses of the shells are estimated. For U Cam 
it is found to be 6x10 ''Mq, and for DR Ser 9xlO"^M0. 
Both estimates are found to be consistent with previous re- 
sults. 

- The ages of the detached shells around U Cam and DR Ser 
are < 700 yr and < 1400 yr, respectively. The measured 
widths of the shells imply formation time-scales of a few 
hundred years. This is consistent with the scenario of de- 
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tached shells forming as an effect of thermal pulses and sub- 
sequent wind-interaction. 
- The CO shell around U Cam, as measured by previous inter- 
ferometric observations, appears smaller than the dust shell 
imaged in this paper, indicating a possible separation of the 
dust and gas since the formation of the shell. 
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Appendix A: The PolCor instrument 

The main reason for the construction of PolCor was the need 
for an instrument that can measure faint scattered light close to 
bright stars. The instrument is optimized for both scattered light 
from circumstellar dust particles as well as resonance line scat- 
tering from circumstellar gas in the following ways. To increase 
the contrast between the PSF wings and polarized scattered light 
from the dust, the instrument includes a polarizing mode. To fur- 
ther bring down the surface brightness of the wings of the stellar 
PSF (Point Spread Function) and avoid saturation of the central 
star, a coronographic optical design was chosen. In addition, to 
cancel out the diffraction cross of the PSF, a Lyot stop blocks the 
image of the support blades of the secondary mirror. In order to 
optimize the contrast ratio in the detection of resonance line scat- 
tering, the instrument is equipped with ultra-narrow band optical 
filters. Finally, to spatially resolve structures in the circumstel- 
lar environments, the instrument uses Lucky imaging (Sects 13.21 
and lA.2b . which considerably improves the sharpness of the im- 
ages compared to the seeing limited case. 

A. 1 . Technical details 

The PolCor instrument is briefly described in the following 
points: 

1. The EMCCD (Electron Multiplying CCD) camera (Andor 
IXON) uses a thinned 512x512 CCD array with 16 pix- 
els, giving a full field-of-view of 1 arcminute. For low light 
levels it can be used in a photon counting mode, which is in 
principle noiseless. In practice this mode is limited by clock 
induced pulses, occurring typically once (per pixel) for 200 
readouts. The sky emission is too bright for photon count- 
ing with broad band filters, and the normal mode of oper- 
ation is the EM mode. Our laboratory measurements show 
that the EM mode is linear over a very wide brightness range 
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(5 orders of magnitude!). In practice, all frames are stored 
so the mode of operation is a post-processing decision. The 
fastest full-frame readout rate is 33 Hz. Very high time res- 
olution can be achieved by limiting the readout area of the 
chip, making speckle interferometry possible. The quantum 
efficiency is not particularly high (around 30%) in the UV, 
but this is probably because the anti-reflection coating is de- 
signed for longer wavelengths. The efficiency of the anti- 
reflection coating is high enough in the red spectral region 
that no interference fringes (caused by the OH sky emission) 
have been seen in the observations. 

2. A high-quality polarizer (Meadowlark Optics, DP-050-VIS) 
is placed in the converging beam from the telescope and can 
be rapidly turned to the four different positions: 0, 45, 90 
and 135 degrees. At each position, typically 30 images with 
0.1s integration time are taken. One unit measurement cy- 
cle (which also includes a dark measurement with a closed 
shutter) takes approximately 20 s. It is repeated 100 times 
in order to compensate for changing sky conditions (seeing 
and transmission). The EMCCD has a frame transfer read- 
out, which means that no time is lost due to readout. As a 
consequence, the overhead is limited to the time for turning 
the polarizer and the resulting overhead is 20% of the on tar- 
get time. 

3. The coronographic masks, positioned in the focal plane of 
the telescope, consist of neutral density (ND) disks with 
three different sizes corresponding to 1'.'5, 3", and 6" on the 
NOT. For each size, three different ND:s are available: ND - 
2, 3.5 and 5, corresponding to damping factors of 100, 3300, 
and 10^. The main reason for not using opaque disks is the 
need for an accurate centering of the star This can be conve- 
niently achieved when the star can be seen through the mask. 

4. The re-imaging optics are based on mirrors which provide a 
diffraction limited performance. To achieve the same image 
quality with lens optics, one would have to use several ele- 
ments in both the collimator and the camera, making it hard 
to avoid disturbing ghost images due to multiple reflections. 
The reflectivity of the coatings (CVI Laser Corp.) of the four 
mirrors (two off-axis paraboloids and two flat folding mir- 
rors) is better than 98.5% per surface over the whole sensi- 
tivity range of the detector, and thus the level of scattered 
light is kept low. The re-imaging optics have a 1 : 1 magni- 
fication (which is equivalent to a pixel scale of 0'.'12 at the 
NOT). Two Barlow lenses are available for 2x and 3x mag- 
nification (intended for speckle interferometry). 

5. To avoid the diffraction cross from the secondary mirror sup- 
port of the telescope, the Lyot stop blocks the re-imaged 
cross (with wider bars) as well as the secondary mirror 
(slightly oversized). Due to the Alt-Az construction of the 
NOT, the field de-rotator causes the diffraction cross to ro- 
tate and to compensate for that, the Lyot cross has a com- 
puter controlled de-rotator 

6. Standard broad-band ^/feri (Bessel U, B, V, R and I) are 
available as well as ultra-narrow band (x; 1 A) filters for the 
Call, Nal and KI resonance lines. For each of the reso- 
nance line filters, double-peaked reference filters are avail- 
able. The filter holder flips quickly (x; 1 s) between two po- 
sitions, which allows for multiple filter exchanges during 
a measurement. This allows for accurate observations of 
line/continuum ratios also during less good photometric sky 
conditions. 



A.2. PolCordata reduction and performance 

As a first step in the data reduction, the image motion, defined 
by the centroid of the light distribution of a stellar image, is de- 
termined. Figure [ATI illustrates the image motion for a period 
of good seeing conditions. The images show two field stars. The 
star to the left is brighter than the star to the right. The eight 
images in Fig. lA. H are taken at an interval of 0.1 s and the plus 
sign denotes the position of the brighter star for the average im- 
age of 1500 exposers. Figure [AS] shows the equivalent observa- 
tions during worse conditions. Even at good seeing conditions, 
the image motion is quite noticeable. In Fig. lA.ll the speckles 
overlap, while in the period of poor seeing (Fig. IA.2) they are 
more or less spread out. The image motion is typically a small 
fraction of an arcsecond on short time-scales. The image motion 
is further illustrated in Fig. \A.3\ (lower panel) where the relative 
center positions for 1500 frames (150 s) observed during good 
conditions are shown. A frame rate of 10 s"' suffices to resolve 
the image motion {upper panel). 

In order to quantify the spread-out due to the image mo- 
tion, we calculate the sharpness of each image. The sharpness 
is measured as the percentage of light that enters a box with 
sides - (y.'55 centered on the brightness peak. It varies on short 
time scales, in particular during periods of poor seeing (Figs lA.2l 
and lA.4b . 

By correcting for the image motion, the resulting seeing is 
usually improved by 20-30%. By calculating the light concen- 
tration for each image and only including a fraction of them in 
the shift-and-add reduction step, the final image can be much 
sharper than that for traditional long integration imaging. It 
should also be noted that tracking errors of the telescope are 
easily compensated for in this shift-and-add reduction step. As a 
trade-off between sharpness and depth an acceptance level needs 
to be chosen, and shift-and-add procedure is then performed only 
for the accepted frames. In Fig. lA.5l the effects of only shift-and- 
add and selecting 15% of the frames are shown. 

The main purpose of PolCor is to measure faint scattered 
radiation close to bright stars and the level of diffracted and 
scattered light of the telescope/instrument combination should 
be as low as possible. In Fig. lA.6l the PSF is shown for both the 
cases with and without an occulting disc. 
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Fig. A.2. Same as Fig. lA.ll but observed under poorer seeing conditions. Each image shows a pattern of more or less scattered 
speckles. The sharpness (measured as the percentage of the light that enters a box with sides = 0.55") is indicated in the upper right 
corner of each image. 
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Fig.A.4. The sharpness of an image varies marginally during 
good seeing conditions (upper curve, corresponding to the exam- 
ple shown in Fig. lA. Il l, and rapidly during poor seeing conditions 
(lower curve, corresponding to the example shown in Fig. IA.2b . 
The frame rate is 10 s"' and the sharpness is measured as the 
percentage of the light that enters a box with sides - 0.55". 
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Fig. A.3. Upper panel: The image motion during 3 s. The time 
line is color coded to make it easier to follow. The used frame 
rate (10 s"') clearly suffices to resolve the image motion. Lower 
panel: The center positions for a star during 150 s (1500 frames). 
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Fig. A.5. The seeing in images attained by simply averaging all 
the 1500 frames in the two examples observed under good {up- 
per panel) and poor conditions {lower panel) is represented by 
the red curves. The FWHM is 0.7" for the upper and 1.1" for the 
lower curve. The blue curves show the results when only shifting 
and co-adding the frames. The black curves represent the seeing 
when only the sharpest 15% of the frames are used for the co- 
added image. The FWHM is 0.4" for the upper case and 0.7" 
for the lower, clearly demonstrating the improvement in spatial 
resolution. 
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Fig.A.6. Upper panel: The radial cut of the PSF for PolCor 
without any occulting disc. The surface brightness is given as 
the fraction of the total intensity of the star per square arc sec- 
ond. Lower panel: The radial cut of the PSF is shown for the 
case when an occulting disc with a diameter of 6" and attenua- 
tion of a factor 100 is used. 



